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STABLE BLOWUP FOR WAVE EQUATIONS IN ODD 
SPACE DIMENSIONS 


ROLAND DONNINGER AND BIRGIT SCHORKHUBER 

Abstract. We consider semilinear wave equations with focusing power 
nonlinearities in odd space dimensions d > 5. We prove that for every 

. , o d+1 d — 1 

p > there exists an open set of radial initial data in if 2 x H 2 
such that the corresponding solution exists in a backward lightcone and 
approaches the ODE blowup profile. The result covers the entire range 
of energy supercritical nonlinearities and extends our previous work for 
the three-dimensional radial wave equation to higher space dimensions. 


1. Introduction 


We consider the initial value problem for the focusing nonlinear wave equa¬ 
tion 


dfu — An = \u\ p 1 n, 
u\ t =o = u 0 , d t u\ t =o = ui, 


( 1 . 1 ) 


for (t, x) £ I x M d , d = 2k + 1, k > 2 and I an interval, where 0 £ I. 
Eq. m is conformally invariant for p = and we restrict ourselves to 

the superconformal case 


p > 


d + 3 

d- 1' 


( 1 . 2 ) 


The above equation enjoys scaling invariance in the sense that if u solves 
Eq. G3D then another solution can be obtained by setting u\(t,x) := 

2 

A p- 1 u(t/X, x/X) for A > 0. The conserved energy is given by 

E(u(t,-),d t u(t,-)) = i||(n(t,-),^n(t,-))||^i xL2(]Rd) - ^IK*, ■)II5 +i ( r«* ) 


and it is invariant under the above scaling for p = which defines the 
energy critical case. In general, the scaling invariant Sobolev spaces are 
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H Sp x H Sp 1 (M d ), where the index s p = | is usually referred to as the 

critical regularity. 


1.1. Basic well-posedness theory and explicit blowup solutions. 

One is usually interested in (strong) solutions of Eq. (11.11) that satisfy the 
equation in integral form by using Duhamel’s principle, see for example [45] . 
In this sense, Eq. (11.11) is locally well-posed in H Sp x H Sp ~ 1 ( R rf ) for d > 5 
and p > jrfj, given that the nonlinearity is sufficiently regular, cf. Lindblad 
and Sogge j32j. Moreover, solutions that correspond to sufficiently small 
initial data can be extended globally in time. We also note that local well- 
posedness in H s x H s ~ 1 (R d ) for s > § and smooth nonlinearities is classical 
[35] . However, global well-posedness does not hold in general. A convexity 
argument by Levine m shows that initial data with negative energy (and 
finite L 2 —norm) lead to blowup in finite time, cf. also [23] for generalizations. 
Explicit examples for singularity formation can be obtained by considering 
the so called ODE blowup solution 


u T (t,x) = c p (T - t) p- 1 , 


Cr) • — 


2(P+1) 


1 

P-1 


(1.3) 


which is independent of the space dimension and solves the ordinary differ¬ 
ential equation utt = \u\ p ^ 1 u for p > 1. By finite speed of propagation one 
can use ut to construct compactly supported smooth initial data such that 
the solution blows up as t —>• T. 

In one space dimension the ODE blowup mechanism is universal, cf. the 
fundamental work by Merle and Zaag [mi, [Mi, m, E 2 ] and the references 
therein. In higher dimensions, the situation is more complex. Depending on 
d and p many other explicit examples for singular solutions were found in 
the past years, including the celebrated work of Krieger, Schlag and Tataru 
[29j on type II blowup solutions for the energy critical equation in three 
space dimensions, see below. For d = 3, p = 3 and p > 7 an odd integer, 
it was proved by Bizon, Breitenlohner, Maison and Wasserman 0 , m that 
Eq. (11.11) admits infinitely many radial self-similar blowup solutions of the 

form ( T—t ) p- 1 n € No, with ut corresponding to the groundstate, 

i.e., /o = c p . Another blowup mechanism for Eq. (ED, which only exists 
for d > 11 and a range of supercritical nonlinearities p > p(d) > was 
recently established by Collot [6j, see below. 

Most of these explicit solutions have unstable directions, i.e., they are un¬ 
stable under generic small perturbations and are not supposed to describe 
the ’typical’ blowup behavior for solutions of Eq. ED, see for example [26] . 
On the other hand, numerical experiments by Bizon, Chmaj and Tabor [2] 
for the three-dimensional equation show that the behavior of generic radial 
blowup solutions can be characterized in terms of the ODE blowup solution 
locally around the blowup point. 

The stability of ut in three space dimensions was established in our previous 
works Ena, m for radial perturbations and all p > 1. Recently, we could 
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extend this to the general case (without symmetry) {TTJ for p > 3. For 
subconformal nonlinearities the dynamics around ut were also investigated 
by Merle and Zaag m in the non-radial setting and in arbitrary space 
dimensions. 

In view of the findings in |6j for supercritical radial wave equations in high 
space dimensions, we extend our previous results and establish the stability 
of the ODE blowup solution in arbitrary odd space dimensions. Although 
for d = 3 we were able to drop the symmetry assumption, it is an open 
question how this can be accomplished for d > 5, see the discussion below. 
We therefore restrict ourselves to the radial case and study solutions that 
blow up at the origin (which is the most interesting case). 

We note that this work is not a mere technical generalization of [10] , H2]. it 
can rather be viewed as a systematization and refinement of our approach 
that has also been applied (with slight modifications) to establish stable self¬ 
similar blowup for equivariant wave maps P3], [JO] and Yang-Mills fields [8] 
in supercritical dimensions. 

1.2. Radial solutions in lightcones. In the following we use the abbre¬ 
viation u[t] := ( u(t , -),dtu(t, •)). We are interested in the behavior of radial 
solutions of Eq. (11.111 in backward lightcones 

r T (M d ) := {(t,®) € [0, T) x R d : |s| < T-t}, 

with vertex (T, 0) for T > 0. Consequently, a suitable concept of (strong) 
solutions in lightcones is required. This can be obtained for example by 
combining the classical Duhamel formula on R d with suitable cut-off tech¬ 
niques, see [21] • Here, we pursue another approach which is based on the 
formulation of Eq. (11.11) in self-similar coordinates 

t := — log(T-f) TlogT, £ 

To motivate the following, let u € C' 0O (Tt’) be a radial solution of Eq. (jl.ll) . 
By setting 

^(-log(T-f)+logT,^) := (T -t)~u(t,x), 

ft(~^g(T - t) YlogT, := (T - t)~ +i d t u(t,x), 
we obtain a smooth solution of the first order system 

drift = ft - £ • S/ft - ^zrft, 

d T ft = Aft - £ • s/ft - 2±iV’ 2 T + ft\ ftr 1 . 1-5 

By setting T 7 (r) := (ft (r, •), ft (r, •)) this can be written as 
9 t ^ t (t) = /: 0 ^ t (t) + F(^ t (t)), 

where represents the linear part of the right hand side of Eq. (11.51) . To 
formulate the following statement we define for A: € No 

H^ d (M d ) :={ue H k (M d ) : u is radial}. 
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d+1 d—1 

Proposition 1.1. Let f) := H rad x H rad (B rf ). There is a dense domain 
T>(Co) C S) such that the operator Cq : D(£o) C S) -A S) is closable and its 
closure generates a strongly continuous one parameter semigroup {5o(r) : 
r > 0} of bounded operators on S). 

This is an immediate consequence of the results proved in Section 14.2.11 By 
Duhamel’s principle we can now formulate the above equation as an abstract 
integral equation 

T t (t) = S 0 (t)T t ( 0) + r S 0 (r - tOF^V))^'. (1-6) 

Jo 

We take this as a defining equation for our notion of strong lightcone solu¬ 
tions. 

Definition 1.2. We say that u: T^R^) -A R is a radial (H “a - -)solution 
of Eq. (11.11) if the corresponding 4 ’ T belongs to C([0, oo), Sj) and satisfies 
Eq. (11.61) for all r > 0. 

d+l d—1 

Definition 1.3. Let ( 1 / 0 , 111 ) € JJ rad x ff rad (R d ). We say that T > 0 
belongs to T{uq,u\) C (0,oo) if there exists a solution u: rr(R d ) A R to 
Eq. m with it[0] = (rto,ni)| B d . Set 

T (u 0 , ui ) ■= sup{T(u 0 ,ui) U {0}}. 

If T( UOiU1 ) < oo, we call T/ U0>U1 \ the blowup time at the origin. 

1.3. The main result. We prove the stability of the ODE blowup solution 
in the following sense. 

Theorem 1.4. For d = 2k + 1, k € N, k > 2, fix p > and Tq > 0. 
There are constants M, S > 0 such that if uq , u\ are radial functions with 

\\(u 0 , ui ) - utM\ h ^± xH ^ {Ko+s) < -M (1-7) 

the following statements hold: 

i) The blowup time at the origin T := T^ uo U1 ) is contained in the interval 
[T 0 -S,T 0 + S\. 

ii) The solution u : r^R^) —> R satisfies 

(T-t)F^\\u{t, ■) -u T {t,-)\\ L oo {M d_ t) < ( T-t ) Mp , 

(T - t)- s r\\u(t, •) - u T (t, •)llL 2 ( B d_ t ) < (T - t )>•*, (1.8) 

{T - t)~ s r +1 \\u[t] - u T [t]\\ tilxL2{B d T _ t) < (: T-tY p , 

for s p = |- Up = minj-^j, 1} — e and e > 0 small. Furthermore, 

forj = 2,...^, 
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Remark 1. The normalizing factors in Eq. (11.81) and Eq. (|1.9j) appear natu¬ 
rally and reflect the behavior of ut in the respective norms. Since the ODE 
blowup solution has a trivial spatial profile, it vanishes identically in higher 
order homogeneous Sobolev norms, which yields Eq. (ll.Djl . The e-loss in 
the convergence rates is due to the application of abstract arguments from 
semigroup theory. 

Remark 2. Our approach is perturbative, i.e., we construct solutions of the 

d+1 

form u = ut + This and the embedding H 2 ^ L°° guarantee that the 
nonlinearity is smooth for all p > 1 provided that the perturbation is small 
enough. In particular, Theorem 11.41 can be extended to all p > 1 without 
modifications. We are therefore able to construct solutions of Eq. (11.11) (at 
least in a backward lightcone) for nonlinearities that are not covered by the 
standard local well-posedness theory. 

If we restrict ourselves to Sobolev spaces of integer order, then the regularity 
required in Theorem 1 1.41 is optimal for p > 5 by local well-posedness ([s p ] = 
^33). However, for p < 5 this might be improved and we show this explicitly 
for p = 3 in Theorem [L5] below. In this case, [S 3 ] = and If ^-solutions 
can be defined in a similar manner as above. 

Theorem 1.5 (Improvement of the toplogy). Let p = 3. Then the first 
statement of Theorem \l-4\ holds (mutatis mutandis) for radial initial data 
satisfying Eq. (|1.7I) in H^ Sp ^ x , where s p = | — 1. The corresponding 

solution u : Tt^M^) —>• R satisfies the estimates 

(T-t)~ Sp \\u(t,-) - u T {t,-)\\ L2(M d_ t) < (T-f)5- £ , 

(T -t)- Sp+] \\u[t] ~ U T [t}\\H jx Hj-l (M d_ t) < (T-f)5- £ , 

for j = 1,..., [s p ]. 

We note that there are possibly other situations where the topology can be 
optimized and it will become clear in Section 11.51 how this could be realized 
within our framework. However, we do not pursue this here. 

1.4. Related results. Blowup for the wave equation with (sub)conformal 
focusing nonlinearities 1 < p < 4^, d > 2 , was considered in the seminal 
work of Merle and Zaag [33], [35] , cf. also Antonini and Merle p]. They were 
able to prove that all blowup solutions diverge with the self-similar rate in 
the backward lightcone of the blowup point (at the energy level). They also 
extended their previous analysis for the one-dimensional wave equation to d 
dimensions in the subconformal case to characterize the behavior of radial 
solutions provided that the blowup occurs outside the origin [38]. In m, 
[42 ] they studied the dynamics around ut (without symmetry assumptions) 
and investigated properties of the blowup surface. 

In the superconformal regime, much less is known concerning the behavior 
of generic solutions. However, for energy subcritical nonlinearities, Killip, 
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Stovall and Visan [23] as well as Hamza and Zaag [19] were able to derive 
upper bounds for the blowup rate. 

In the energy critical case, ut is the unique self-similar solution and it can be 
used to construct blowup solutions that diverge in the scale invariant norm 
H 1 x L 2 (R' 3 ), cf. for example [30]. This behavior is referred to as type I and 
contrasted by type II blowup, where solutions stay bounded in the critical 
norm. First examples of type II solutions were obtained by Krieger, Schlag 
and Tataru m, m for the radial equation in three dimensions using the (up 
to scaling) unique solution of the corresponding elliptic problem. We also 
refer to [9] for solutions that blow up in infinite time. A detailed description 
of all possible type II blowup dynamics was provided in the celebrated work 
of Duyckaerts, Kenig and Merle, cf. |14| as well as m for further references. 
For d = 4, smooth type II blowup solutions were constructed by Hillairet 
and Raphael [20], and we also refer to [21] for a recent result by Jendrej in 
five dimensions. 

To the best knowledge of the authors, all currently available results for 
supercritical nonlinearities p > are either conditional or consider per¬ 
turbations around certain special solutions. Type II blowup behavior for 
radial solutions was excluded by Duyckaerts, Kenig and Merle [T5] for d = 3 
and by Dodson and Lawrie [7] for d = 5. We also refer to similar results for 
the defocusing case obtained in [23], [25] or [5]. These works show that if 
the critical norm stays bounded up to the maximal time of existence, then 
the solution is global in time. We also mention a recent work by Krieger 
and Schlag [28], where smooth global solutions are constructed that have in¬ 
finite critical norm (but are bounded in all higher norms). Recently, based 
on the pioneering work of Merle, Raphael and Rodnianski [33] for the en¬ 
ergy supercritical nonlinear Schrodinger equation, a new blowup mechanism 
was described by Collot [6] for the radial wave equation in d > 11 and 
p > 1 + ——| , . There, solutions blow up via concentration of the soli- 
ton profile, which is somewhat reminiscent of the type II behavior in the 
energy critical case. However, the solutions diverge in the critical norm and 
this blowup mechanism could therefore be referred to as type lib in order 
to avoid confusion. 

1.5. Strategy of the proof. We consider the radial equation 

dtu - ^rn:d r (r d ~ l d r u) = |u| p-1 u, (1.10) 

for u = u(t , r ) and initial data (uq, u\) = ut 0 [0] + (/, g ), where (/, g) are free 
radial functions. We study the initial value problem in a backward lightcone 
{(t, r) : t E [0, T),r € [0, T — f]}, where T > 0 is a parameter that will be 
fixed in the final step of the proof. We introduce rescaled variables and 
rewrite Eq. (11.101) as a first order system in (radial) similarity coordinates 

— log(T ~ t) + log T , 


r := 
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which yields the abstract evolution problem 

dr^ = Lo'h + F(T), 

where T = (Vh,^)- Here, Lo represents the radial wave operator and 
F('k) = (0, |-01 | p—1, 0i)- The backward lightcone now corresponds to 

{(r,p) : t € [0, oo),p € [0,1]}. 

Note that the parameter T does not appear in the equation itself, but it 
shows up in the initial data. In this formulation, the ODE blowup solution 
corresponds to the static solution c p . The ansatz T = c p + <h yields 

d T <I> = L<h + N(<f>), (1.11) 


with L := Lo + L/ representing the linearized part of the equation and N 
denoting the nonlinear remainder. Eq. (11.111) is investigated as an abstract 
ODE on a Hilbert space with norm 


u 


IMI ■ 1)11 (is d ) t II^G 


\ 2 H k-l 


for u = (ui,ti 2 ). For Theorem 11.41 we choose k = (d + l)/2, whereas for the 
proof of Theorem 1 1.5 1 k = {d— l)/2 is sufficient. The choice of the topology 
is motivated as follows. 

First, we have to derive suitable Lipschitz estimates for N. In the first situ¬ 
ation (Theorem 11.41) . we can exploit the Sobolev embedding H k L°° for 
k > d to infer that the nonlinear remainder is smooth for all p > 1 (given 
that the perturbations are sufficiently small). An application of Moser’s 
inequality then yields the desired result. For p = 3 (Theorem 11.51) . the 
nonlinearity is analytic, hence regularity is not an issue and the Lipschitz 
estimates can be obtained by using Holder’s inequality and Sobolev embed¬ 
ding. 

Furthermore, we need a decay estimate for the time evolution of the linear 
wave equation in similarity coordinates. To see what can be expected, let 
us drop the symmetry assumption for a moment and let u(t,x), x € M d , be 
a generic solution of the free wave equation 


dfu — A u = 0. 


Let ipi denote the rescaled solution in similarity coordinates as in Eq. (hd. 
The scaling behavior of Sobolev norms implies that 

IIVhXD = — t) p_1 ^ 


where we write r = — log(T — t) + logT for brevity. One can easily check 
that \\u(t, ■)|| i jfc( B d ) is bounded. However, without further assumptions 
on the regularity this cannot be improved since one can construct explicit 
solutions that decay arbitrarily slow in H k { Hence, a decay estimate 
for tjji in H k (M d ) can only be obtained if | + k > 0. Unfortunately, in 
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backward lightcones such homogeneous quantities are only seminorms and 
we have to work instead with 

k 

3=0 

At first glance, the lower order terms seem to spoil the decay estimate. 
However, this can be overcome by considering equivalent norms. In the 
radial case, the construction is based on the reduction of the d— dimensional 
radial wave equation to the one-dimensional case (or simply to a lower 
dimensional equation, depending on the required level of regularity). For 

our purpose we set Du(t,r) := (~<9 r )^~ (r d ~ 2 u(t,r)) and observe that 
(d-i )/2 

(T — t)p- 1 2 J ||u(i, | ■ l)llijj(Bd_ t ) 

3=0 

— ( T — t)p- 1 2 \\Du(t, ')lli/i(o,T_q' 

Since Du solves the one-dimensional equation, \\Du(t, Olljjim T—t) bounded. 
This equivalence is crucial in the proof of Theorem 11.51 It is also obvious 
that we have decay only if p < 5. For Theorem ll.4l this is not sufficient (also 
because the above quantity only provides derivatives). We could work 
instead with \\Du(t, -)||p- 2 ( 0 t)’ but this is only a seminorm (the radial 
derivative of Du does not vanish at the origin). Adding the energy part 
solves this problem, but spoils again the decay estimate for p > 5. Hence, 
we consider 

■= (T-i)i^ +1 ||Du(t,0||^ (OiT _ t) 

+ (T-t)^ +1 \\Dut(t,-)\\% HOiT _ t) 

+ (T - t) Wi | (Du) r (t, T — t) + ( Du)t{t , T-t) I 2 

and show that II • || d+ 1 d- 1 ~ || • ||n- Observe that solutions of the 

" "ifTxU-2“(B d ) " " 

one-dimensional equation wtt — w rr = 0 satisfy 

jt\w r (t,T - t) + Wt(t, T — t)\ 2 = 0. 

In view of this, the desired decay follows for all p > 1. We note that 
for d = 3 such an equivalent norm also exists in the non-radial context, 
cf. m- However, it is not clear how this can be generalized to arbitrary 
space dimensions. 

To prove the results of Section 11.31 we proceed as in (10], m and use the 
theory of strongly continuous one-parameter semigroups to address the lin¬ 
earized equation. Since the operator L is a highly non-selfadjoint object, 
semigroup theory can deploy its full strength and enables us to treat the 
problem on a very abstract level. 
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With the above considerations and suitable equivalent norms it is 
easy to show that Lo is the generator of a semigroup which satisfies 
a suitable decay estimate. Since L/ is bounded, well-posedness of 
the linearized problem and the existence of a strongly continuous 
semigroup (S(r)) T >o generated by L follow immediately. 

To deduce suitable growth estimates for the semigroup we analyze 
the spectrum of the generator. Compactness of the perturbation re¬ 
duces matters to the investigation of the eigenvalue problem, which 
can be solved explicitly in terms of hypergeometric functions. We 
show that the spectrum of L is contained in a left half plane ex¬ 
cept for the point A = 1, which is an eigenvalue with eigenfunction 
g. The existence of this unstable eigenvalue is a consequence of the 
time translation symmetry of the problem and we define a spectral 
projection P to analyze the behavior of solutions on the stable sub¬ 
space. Note that we have to verify that rgP = (g), since we are 
dealing with a non-selfadjoint problem. In contrast to m , m, 
where we used resolvent estimates and the Gearhardt-Priiss Theo¬ 
rem to deduce growth bounds for (1 — P)S(r), we employ a much 
simpler argument here and exploit the compactness of the pertur¬ 
bation directly. This is a substantial simplification relying only on 
standard results from semigroup theory. As a result we obtain that 

||(1 — P)S(t)|| < e -MpT , PS(t)u = e T u, where fi p > 0. 

We rewrite Eq. ( 11,111) in Duhamel form, 

$(t) = S(t)U(v,T) + [ S(T-r , )N($(r'))dT / , (1.12) 

Jo 

where U(v,T) with v = (/, g) gives the original initial data. The 
main ingredients for the nonlinear theory are the above estimates 
for the semigroup and Lipschitz estimates for the nonlinearity of the 
form 

||N(u)-N(v)|| < (||u||+ ||v||)||u-v||. 

The rest of the proof is purely abstract. 

We add a correction term to Eq. (11.121) in order to suppress the 
unstable behavior of S(r). An application of the Banach fixed point 
theorem shows the existence of a unique solution to the modified 
equation 

$(r) = S(r)U(v,T)+ [ S(r - tON^t'))^' 

Jo 

-e T C($,U (v,T)), 

given that v is small and T is close to To. Furthermore, the solution 
decays to zero with the linear decay rate. 

In the final step, we show that for every small v there exists a T v 
close to Tq such that C(<3?, U(v, T v )) = 0. We exploit the fact that 
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C($,u) G (g) and apply the Brouwer fixed point theorem as in m- 
This is a substantial simplification compared to m,n where dif¬ 
ferentiability of several quantities was required. Transforming back 
to original coordinates yields the result. 

2. Notation 

Throughout the paper we assume that d = 2k + 1, k G N, k > 2 is fixed. 
We write N for the natural numbers {1, 2, 3,... } and set No := {0} U N. 
As mentioned above, denotes the open ball in M. d centered at zero with 
radius R > 0. If R = 1, we simply write M d . The notation a < b means 
a < Cb for an absolute constant C > 0 and we also write a ~ b if a < b and 

b< a - ^ 

For a function x i-a g(x), we use the notation g^ n \x) = d for derivatives 
of order n G N. For n = 1,2 we also write g'(x) and g"(x), respectively. For 
a function (x,y) i-A f(x,y) partial derivatives of order n will be denoted by 
dxf{x,y) = tfgnf{x, y) = dff(x,y). For Q C R d a domain, H m (n), m G N 0 , 
denotes the standard Sobolev space with norm 

IMI 2 h™(Q) := 

a:\oi\<m 

where a = (a\,... ,ad) G Ng is a multi-index, i.e., d a = and 

| ck | — cx\ + • • • + a d . 

The set of bounded linear operators on a Hilbert space R is denoted by B(R). 
For a closed linear operator L we write cr(L) and cr p (L) for the spectrum 
and point spectrum, respectively. Furthermore, we set Rl(A) := (A — L) -1 
for A ^ <r(L). 

3. The radial wave equation in similarity coordinates 

We restrict ourselves to radial solutions of Eq. m and write u(t, x) = 
u(t,r ), where r = |x|, by slight abuse of notation. We introduce the radial 
Laplace operator on M d , 

A r u(t, r) := d^u(t, r ) + d r u(t , r) 

and study the equation 

dfu(t, r) — A r u(t, r ) = | u(t, r)\ p ~ 1 u(t, r ). (3.1) 

The initial data at t = 0 are assumed to be of the form 

«[0] =u To [0] + (f,g), (3.2) 

where To > 0 is fixed and (/, g) can be chosen freely. At the origin we 
impose the natural boundary condition d r u(t, 0) = 0 for all t > 0. We define 
rescaled variables 

Ui{t,r) := (T - t)p-iu(t, r), U 2 {t,r) := (T - t) 1+ ^d t u(t,r), 
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where T > 0. This yields the first order system 

\(T- t)Arlh - (T - t)- 1 (0U 2 - pt^Ui) 


d t U i 
d t U 2 


with data 


— _ p+ 1 

Ui(0,r) = (^-) — Cp + T~f(r), U 2 (0,r) = ^jC p + g ( r ). 

and boundary conditions 

d r Ui(t, 0) = d r U 2 (t , 0) = 0, 

for all t > 0. In the rescaled variables the blow-up solution ut is static and 


P+1 


T = - log (T - t) + log T. 


corresponds to c p := (c p , -c p ). We introduce similarity variables 

r 

P ~ T-t ’ 

Derivatives transform according to 

+ M>), d r = e —dp. 

Setting ^-(r, p) := Uj(T — Te~ T ,Te~ T p) for j = 1,2, we obtain the system 


d T ^i 

2 


V>2 - pdp^l - ^+ T "01 

A p Vh - pdpi >2 - |rT'i/’2 + |Vh| p_1 Vh 


(3.3) 


with the boundary conditions <++,■ (r, 0) = 0 and initial data 


(3.4) 


V’l(0,p) = (J+-iCp + Tr-i/(Tp), 

m P+1 0 P+1 

p) = (^)r- 1 ? - T C p + Tr-i 5 (Tp). 

We restrict the problem to the backward lightcone of (T, 0), i.e., we study 
Eq. (j3.3f) for p € [0,1] and r > 0. 

3.1. Perturbations around the ODE blow up solution. Inserting the 
ansatz 


Vh 

V>2 


<£>2 


+ Cr 


into Eq. (13.311 we obtain 

<9 t </21 
d r <P2 


ip 2 - pdpipi - ^(pi 
PA 
P~ 


yAp^i - pd p ip 2 - +pcp Vi + A^i) 

for p G [0,1] and r > 0 where 

N{tp 1 ) = | c p + V9i| p_1 (c p + (pi) - (? p - p(? p - l p 1 . 

The initial data are given by 

¥>i(0 ,/») = (^)“c p + T“/(Tp)-Cp 

P+1 P+1 

^ 2 ( 0 , p) = (^)T rT ^ T Cp + TWT 5 r(T / 9) - ^-Cp. 
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4. Proof of Theorem 11.41 
Throughout this section, p > is a fixed real number. 

4.1. Functional setting. We consider radial functions u defined on K 
i.e„ «(f) = u(|£|), £ € K- In order to avoid confusion owing to identification 
of u and u, we define 

H r r n (Mf { ) :={u : (0, R) —>• C such that u is m — times 

weakly differentiable and ||u(| • |)|| H m {M d R ) < °°}> 

for rn € No. The density of C°°(M R ) in H m (M R ) implies the density of 
C™[0,R] := {u € C°°[0,R] : i/ 2fc+1 )(0) = 0,Jfe G N 0 } 
in H™(M d R ). For the rest of the paper we set 

d + 1 

md := ^~ 

and introduce the Hilbert space 

n := H™ d x H™ d ~ l (M d ) 

with norm 

M IMI ' 1)11 (IB d ) 11^2(1 ' 1)11 — 1 (IJrf)’ 

for u = (ui,u 2 ) T - 

4.1.1. Equivalent norms onR. We define a norm on TL which is ’tailor-made’ 
for the investigation of the linearized time evolution of the perturbation. 
First, we need the following auxiliary result. 

Lemma 4.1. Let u € H™ d ( B rf ). Then 

m d 

ll u (l • \)\\ 2 H™d.(B d ) — ll n lli 2 (0,l) + Ys IK')" lv,( ' n ' > lli 2 (0,l)' 

n= 1 

Furthermore, for all u € H r f ld ~ 1 (f B d ), 

m d -1 

IWI ■ DllWw == E ll(')”“ ( ” , lli»(0,D- 

n =0 

The proof is given in Appendix [Bj To proceed, we define 
D d u{p) := (p -1 ^) 2 {p d ~ 2 u{p)). 

Note that 

m d -2 

D d u(p) = Y a n p n+1 u^ n \p) = aipu(p) 4-f V)( p ), (4.2) 

n=0 
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for constants a n > 0. The kernel of D d consists of functions which are highly 
singular at the origin, 


ker D d 


<P- 3 ), 

(p~ 3 ,p- 5 ,...,p-( d - 2) ), 


for d = 5, 
for d > 7. 


We also introduce the integral operator 

d-3 f p 

K d u(p) := p 2 ~ d K~^~u{p), ICu(p) := / su(s)ds. 

Jo 

If u G TL, then u\,u 2 are m.^-times, respectively, (rrid — l)-times weakly 
differentiable functions. By Sobolev embedding, U 2 G C md ~ 2 [5, 1] for every 
0 < 5 < 1 and, since rri ( i > |, u\ € C[0, 1] n C md ~ 1 [8 , 1], Hence, the expres¬ 
sions D d Uj, j = 1,2, are defined as sums of weighted classical derivatives on 
(0,1]. Furthermore, kerU^ = {0} on H™ d (H d ) and H™ d ~ 1 (M d ), respectively. 
We infer that D d is invertible on H and 


K d D d u = D d K d u = u. 

Now, consider the sesquilinear form dehned by 

(u|v) D := (D d ui\D d vi) k2{;dl) + (D d u 2 \D d v 2 ) 

+ {[D d u i]'(l) + [D d u 2 }(l))([D dVl }'(l) + [D d v 2 \(l)), 

and set ||u||£> := y (u|u)£>. For u G H, (D d u\)" and ( D d u 2 )' have to be 
interpreted as sums of weighted weak derivatives. The proof of the next 
Proposition is provided in Appendix [Cj 


Lemma 4.2. || • ||£> and || ■ || are equivalent norms on TL. In particular, 

ll u l | 2 £ ll-Dd«l||fra( 0 i i) + Il-Dd« 2 11^1(0,1) ll u ll D ll u ll 2 
for all u G TL. 

By using the results of Lemma [TT] and Lemma f4.21 together with the Sobolev 
embedding H m { 0,1) C' m ~ 1 [0,1] and the density of C^°[0, l] 2 in TL, we 

obtain the next result. 


Corollary 4.3. Let u G TL. Then D d u\ G C 1 [0,1], D d u 2 G C[0,1] and 

D d Uj(0) = 0, j = 1,2. 

4.2. The linearized problem. We exploit the following commutator rela¬ 
tions satisfied by D d and its inverse K d . 

Lemma 4.4. Let u G Hf nd ~ 1 (B d )riC rnd ~ 1 (0, 1) and let A denote the dilation 
operator defined by A.u(p) := —pu'(p). Then 

D d Au(p ) = A D d u(p) + D d u(p ) (4.3) 

for p G (0,1). Furthermore, for u G H™ d (M d ) n C md ( 0,1) we have the 
identity 


D d A p u(p) = [D d u]"(p) 


( 4 . 4 ) 
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for pe ( 0 , 1 ). 

Proof. We note that the assumptions on u imply that the above expressions 
vanish if and only if u = 0. To prove the identities, we proceed by induction. 
A direct calculation shows that for d = 5, D§Au = AD§u + D§u. Assuming 
that Eq. (14.3j) is true for some odd number d > 5, we use the identities 

(•) 2 Au = A[(-) 2 u] + 2(-)V O-'M' = A[(-)-V] - 2(-)-V, 

to obtain 

D d+2 Au{p) = p~ l [D d [(-) 2 Au]]'(p) = p~ l [D d A[(-) 2 u]]'(p) + 2 D d+2 u(p) 

= P~ l [^ D d[(-) 2 u\]'(p) + 3D d+2 u(p ) = A D d+2 u(p) + D d+2 u(p). 
The identity given in Eq. (14.4p is well-known, cf. j!8j . p. 75. □ 

Lemma 4.5. Let w € C' 1 [0,1] D C 2 (0,1) satisfy w{ 0) = 0. Then 

K d Aw = AK d w - K d w, (4.5) 

and 

K d w" = A p K d w. (4.6) 

Proof. For d = 5, Eq. (14.51) follows from integration by parts. Assume that 
it is true for some d > 5 odd. We use the identities 

KAw = AK.w + 2/C, (-)" 2 Au; = A[(-)~ 2 w] - 2{-)~ 2 w 

to infer that 

K d+2 Aw = (-)- 2 K d JCAw = (• )~ 2 K d AICw + 2 K d+2 w 

= (-)^ 2 A K d JCw + K d+2 w = AK d+2 w - I< d+2 w. 

Using integration by parts one can easily check that Eq. (|4.6I) is true for d = 5 
provided that w( 0) = 0. Assume that it holds for some odd number d > 5. 
To clarify notation we write Ajf^u(p) := p 1 ~~ d (p d ~ 1 u'(p))'. A straightforward 
calculation shows that 

[K d+2 w"]{p) = p~ d [ s d ^ 1 [K d w ,l ](s)ds = p~ d [ s d-1 [A^ K d w\(s)d-s 

Jo Jo 

= p~ 1 [K d w}'(p) = (2 - d)p~ d IC^w(p) + p 2 ~ d IC^w(p) 

= A^K d w(p). 

□ 


4.2.1. Well-posedness of the linearized time evolution. We define the oper¬ 
ator (L 0 ,P(L 0 )) by 


Lqu(p):= 


u 2 {p) + Au x {p) - ^j-ui(p) 
A pUi(p) + A u 2 (p) - u 2 {p) 


(4.7) 
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P(L 0 ) :={uewn C 00 ^, l) 2 : D d u 2 G C 2 [ 0,1], 

D d u\ € C 3 [0,1], [-DdUi] w (0) = 0} . 

Using the results of Lemma 14.41 we get that 

D d (L 0 u)j = (A 0 D d u)j (4.8) 

for j = 1,2, where 

. / V ( w 2 (p) - pw[{p) + ^wi(p) \ 

° W P \ w "(p) - P w 2(p) - ^rW 2 (p) J ' 

In view of Lemma 14.21 is now obvious that the regularity properties satisfied 
by functions in X>(Lq) imply that Lou € PL. We note that 0, l] 2 C T?(Lo), 
i.e., Lo is densely defined. 

Lemma 4.6. Let u E T>{ Lo). T/ien 

Re(L 0 u|u) D < pTu 11 u 11 z) • 

Proof. To abbreviate the notation we set Wj := D d Uj, j = 1,2, for u € 
V(L 0 ). By Eq. (14.81) 

[L> d (L 0 u)i]'(p) + [D d (t 0 n) 2 ]{p) = w 2 (p ) - pw"(p) + w"{p) 

~ pw' 2 {p) - i w i(p) + ^2 (p)) ■ 

Evaluation at p = 1 yields 

[A,(L 0 u)i]'(l) + [A,(L 0 u) 2 ](1) = K(l) + MV] • 

Furthermore, 

[L> d (L 0 u)i]"(p) = w 2 (p) + A w"(p) - 2=zw"(p), 

and 

[D d (L 0 u) 2 ]'(p) = w"\p ) + A w' 2 (p) - Z=jw 2 (p), 

for p E (0,1). Note that for functions w E C[0,1] D C' 1 (0,1) integration by 
parts yields 

Re(Aiu|iu) L 2 (oa) = |||u;||| 2(0]1) - ±Hl)| 2 . 

With this we infer that 

Re([ J D,(L 0 u) 1 ]"|[ J D d u 1 ]") i2(01) = Re KK)^^) ~ \ K'(l)| 2 

~ (s + Fl) IKHi 2 (0,1) ’ 

and 

l 2 (o,i) = Re K'K) L 2 (0i1) - 5 K(l )| 2 


Re([L ) d(L 0 u ) 2 ] , | [D d u 2 ]') 
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Using these identities and performing one additional integration by parts we 
obtain 

Re(L 0 u|u) D = -p^|[iW(l) + [A^ 2 ](l )| 2 

“ + p=l) (||[Ai“l] |lz, 2 ( 0jl ) + II [Pd^ll 11x2(0,1)) 

- ^[D^ni) - [D dU 2 ni )\ 2 

□ 


Lemma 4.7. Fet // = 1 — For every f = (/i,/ 2 ) t € C^fO, l] 2 f/iere 
exists a function u G P(Lq) satisfying the equation 

(/i - L 0 )u = f. 

Proof. Trivially, for f = 0 we have u = 0. Assume that f G C£°[0, l] 2 does 
not vanish identically. We set 

F{p) := D d f 2 {p) + D d hip) + p[D d fi}'(p) 

and define functions 

m{p) := [ 1 2 [ F{s')ds'ds, 

Jo t — s Js 

pi 

w 2 {p) ■■= P 2 / F(s)ds - D d hip)- 
1 - P J P 

The properties of F imply that nq G C 00 ^, 1) fl C 3 [0,1], w 2 G C°°(0,1) n 
C 2 [0,1] and the functions satisfy the boundary conditions ivi(0) = w i (0) = 
w 2 ( 0 ) = 0. A direct calculation shows that wi,w 2 solve the system of 
equations 

pw[ip)-w 2 ip) = Ddhip), .. 

(4.9) 

w 2 (p) - w’lip) + pw' 2 ip) = D d f 2 ip). 


We apply K d to Eq. (14.91) and use the results of Lemma 14.51 to obtain 
K d wiip) ~ K d w 2 ip) - AK'dWiip) = flip), 

2K d w 2 (p) - A p K d wiip) ~ AK'dWiip) = hip). 

Upon setting Ujip) := K d Wj for j = 1,2 and defining u := iui,u 2 ) r we 
obtain a solution of the equation ip — Lo)u = f. The properties of the 
functions Wj imply that u G T>{Lo) and the claim follows. □ 

Lemma 4.8. The operator L' : TL —>• TL defined by 


L'u := 



(4.10) 


is compact. 
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Proof. Let (u n ) n eN be a sequence that is uniformly bounded in PL. By 
Lerrirria f4.2l (DdUi tn ) n ^ is uniformly bounded in H 2 ( 0, 1). The compact em¬ 
bedding H 2 ( 0, 1) H\0, 1) implies the existence of a subsequence, again 
denoted by ( DdU \ tn ) ng ^, which is a Cauchy sequence in H 1 { 0,1). The claim 
follows from the fact that 

||L L Um || ^ pdp \\Dd^l,n Ddlli^mW fll (q jl)- 

□ 

In view of Lemma 14.61 Lemma [4.71 and the boundedness of L', we can apply 
the Lumer-Phillips Theorem mi p. 83, Theorem 3.15, together with the 
Bounded Perturbation Theorem EZj, p. 158, to show that the linearized 
time evolution is well-posed. In particular, by the equivalence of || • ||d and 
|| • || we can formulate the following result. 

Proposition 4.9. The operator (Lo,2?(Lo)) is closable and its closure, de¬ 
noted by (Lo,T ) (Lo)) J generates a strongly-continuous one-parameter semi¬ 
group of bounded operators (So(t)) t >o on TL satisfying the growth estimate 

||So(t)u|| < Me~p~ lT \\u\\ 

for all u £ PL, r > 0 and a constant M > 1. Furthermore, the operator 

L:=L 0 + L', V(L)=V(L 0 ), 

is the generator of a strongly-continuous semigroup (S(r)) r >o. 

In order to derive a suitable growth estimate for S(r) we investigate the 
spectrum of the operator (L,T ) (L)). 

4.2.2. Spectral properties of the generator. 

Lemma 4.10. Let A £ <r(L). Then either A = 1 or 

ReA < max{—^j-, —1}. 

Moreover, A = 1 is an eigenvalue and the corresponding eigenspace is spanned 
by the constant function g = ( 1 , 2 ± 1 ). 

Proof. Let A € <r(L). If ReA < — then the assertion is obviously true. So 
assume that ReA > — Then A ^ cr(Lo) by standard semigroup theory. 
The identity (A — L) = (1 — L'Rl 0 )(A — Lo) and the compactness of L' 
imply that A £ a p ( L). In particular, there exists an eigenfunction u £ T>{ L) 
satisfying the eigenvalue equation (A — L)u = 0. The regularity properties 
of functions in PL imply that L acts as a classical differential operator on 
the interval (0, 1). By a straightforward calculation one can check that 
if u satisfies the eigenvalue equation then u\ £ C[0, 1 ] n C md_1 (0, 1 ) is a 
nontrivial solution of the second order ordinary differential equation 

p 2 u"(p) - Apu(p) + 2(A + )P u '(p) 

+ (A + ^rr)(A + ^=j)u(p) -pc^uip) = 0. 


(4.11) 
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Since the coefficients are smooth on (0,1) we infer that u\ E C* 00 (0,1) n 
C[0,1]. We apply D d to the equation and use the results of Lemma [4.41 
where we proved the identity 

D d[pu'(p)\ = p[D d u}'(p) - D d u{p). 

Similarly, one can show that 

D d [p 2 u"(p)\ = p 2 [D d u]''(p ) - 2 p[D d u]'(p) + 2 D d u(p). 

Upon setting w := D d u\ we infer that 

-(iVKW + 2p(A + ^X(p) , x 

p (4 12) 

+ (A + ^-1)(A + ^Mp)-^-Mp) = 0, 

where w E H 2 ( 0,1) D C 1 )!), 1] satisfies the boundary condition iu(0) = 0, 
cf. Corollary 14.31 By substituting p i-A z := p 2 and setting v(z) := w(y/z) 
one obtains the hypergeometric differential equation 

z( 1 — z)v"{z) + [c — (a + b + 1 )z]v'(z) — abv(z) = 0 

with parameters 

a = |(A-2), 6 = 3 (A + ^f), c=\. 

The assumption ReA > _ implies that Re(c— a — b) = 1 — — ReA < 1. 
Let us assume for the moment that Re(c — a — b) is not zero or a negative 
integer. Around p = 1 two linearly independent solutions are given by 
{ui,hi}, where 

vi (z) = 2 -Pi (a, 6; a + 6 + 1 — c; 1 — z) 

v\ (z) = (1 — z) c ~ a ~ b 2F\{c — a, c — 6; 1 + c — a — 6; 1 — z), 

and 2 -Pi denotes the standard hypergeometric function, see e.g. |33]. If 
Re(c — a — b) = — n, for n € No, then one solution is still given by v\ and 

vi(z) = cvi(z) log(l - z) + (1 - z)~ n h(z), 

where c might be zero for n E N and h is analytic around z = 1. In all cases, 
the requirement w E H 2 (^, 1) excludes the solution v\ and we infer that v 
is a multiple of v±. Around p = 0 we have the fundamental system {uo,uo}> 
where 

vq{z) = 2 F 1 (a,b;c-,z), 

vo(z) = z 1 ^ 2 2 Fi{a + 1 — c, b + 1 — c; 2 — c; z). 

Hence, there are constants co, c\ E C such that 

vi = c 0 v o + ClU 0 . 

The condition iu(0) = 0 implies that u(0) = ui(0) = 0 and thus Co must be 
zero. By [13] . 
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Since the gamma function has no zeros, cq can only vanish if either a + 1 — c 
or 6 + 1 — c = 0 is a pole. This is equivalent to 

A = 1 — 2 k or A = —2k — — ^y for k € No- 

The latter condition implies that A < — ^y which is excluded by assump¬ 
tion. The first condition yields that A = 1 — 2 k for some k € No, hence 

A £ {1, —1, —3, ■ • • 

For p > 3, —£ (—1, 0), hence A = 1 is the only possibility. For < 
p < 3, —oo < —=-r < — 1 and in this case either A = 1 or A < —1. This 
proves the first claim. 

A straightforward calculation shows that g = (1, |i|) satisfies the equation 

(1 - L)g = 0. 

Furthermore, it is easy to check that D^gi = a\p, D^g 2 = C 12 P for constants 
«i, «2 > 0. Hence, g € 'D(Lq), which proves that 1 is an eigenvalue. Suppose 
that there is another eigenfunction g € P associated to A = 1. Then 
Ddgi satisfies Eq. (14.121) . With the same arguments as before we infer that 
Dddi(p ) = • 2 -Fi(a + 1 — c, 6 + 1 — c;2 — c; p 2 ), for some d € C. For A = 1, 

a + 1 — c = 0, hence D<igi = dp which implies that g\ = (3gi for some f3 € C. 
The equation (1 — L)g = 0 then shows that 52 = Pd 2 i which proves that the 
eigenspace of A = 1 is spanned by g. □ 

4.2.3. Time evolution for the linearized problem. 

Lemma 4.11. There exists a projection P G B(TL) onto (g) which commutes 
with S(r) and 

S(r)Pf = e r Pf 

for all f € TL and all t > 0. Moreover, 

||(1 - P)S(r)f || < Me~^ T ||(1 - P)f || (4.13) 

for all f (zTl, t > 0, some constant M > 1 and 

p p = min{^y,l} -e 

for some small s > 0. 

Proof. The eigenvalue A = 1 is isolated and we define P € B(fH) by 

P= 2 ^/ RL(A)dA ’ 

where 7 is a positively oriented circle around 1 in the complex plane with 
radius r 7 = cf. E 21 , p. 178, Theorem 6.5. The projection commutes with 
the operator L and its resolvent, see [22] , p. 173, Theorem 6.5, and thus 
with the semigroup. Furthermore, P. = ker P © rgP and the operator L is 
decomposed into parts L|x>(L)nkerP and L|-p( L ) ni . g p, where L|x>( L ) nrg pU = 
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Lu for u £ D(L) HrgP (analogously for L|-p( L )nkerp)- The spectrum of the 
restricted operator is given by 

cr (T|x)(L)nkerp) = ^(L) \ {1}, cr (I J |x)(L)nrgp) = {li¬ 
lt is immediate that (g) C rgP. It remains to show the reverse inclusion. 
We first observe that if dim rgP = oo, then A = 1 would belong to the 
essential spectrum of L [22], p. 239, Theorem 5.28, which is invariant under 
compact perturbations [22], p. 244, Theorem 5.35. However, 1 0 cr(Lo) and 
we infer that P has finite rank. 

Next, we convince ourselves that rgP C 2?(L). Let u E rgP. By density 
of T>( L) in "H, there exists a sequence (u n ) ne ^ 0 C D(L) such that u n —» u. 
The fact that P is bounded yields Pu„ —» u and since PP(L) C V(L) by 
1221, p- 178, Theorem 6.17, (Pu n ) n£ f} 0 C rgP nP(L). By boundedness of 
L|p(L)nrgP we get tha.t LPu n — >• f, for some f E rgPnP(L). The closedness 
of L now implies that LPu„ —>• Lu and u £ T>( L). We infer that 1 — L| rg p 
acts on a finite dimensional Hilbert space and that A = 0 is its only spectral 
point. Hence, it is nilpotent and (1 — L| rg p) fc u = 0 for all u £ rgP and 
some minimal k £ N. If k = 1, then the claim follows. So let us assume that 
k > 2. Then there exists a nontrivial function u £ rgP C 2?(L) such that 
(1 — L| rg p)u £ ker(l — L| rg p) C ker(l — L), i.e., u satisfies the equation 

(1 - L)u = ag 

for some a £ C. A straightforward calculation shows that the first compo¬ 
nent then satisfies 

P 2 u"(p) - A pUl (p) + ^-pui(p) = JTT a - 

Since v,\ £ C[0,1] fl C md ~ 1 ( 0,1) for u £ H the equation can be interpreted 
in a classical sense for p £ (0,1). Smoothness of the coefficients implies that 
ui € C 00 ]!), 1). We apply D^ and set w := Dju\, where w € 0,1) fl 

C°°(0, 1) and w( 0) = 0 by Lemma 14.31 This yields 

-(1 - p 2 )w''(p) + ^p-pw'(p) - ^p-w(p) = g(p), 

with g{p) = a p p for some a p £ C, a p / 0. Recalling the proof of Lemma 
14.101 we know that a fundamental system is given by w'o(p) = p and 

wi{p) = p 2 ) = (1 -p 2 y^h{p), 

where h is continuous on [0,1] and h( 0) ^ 0 as well as lim p ^i h(p) =4 0. For 

o P + 1 

the Wronskian we obtain W(wo,w\) = — (1 — p z ) p_1 . From the variation 
of constants formula and the boundary condition w( 0) = 0 we infer that 

fp 2 _2_ [P 2 

w(p) = cqp — a v p / sh(s)ds + a p (l — p ) p~ 1 h{p) / s (l — s )p- 1 ds 

J po Jo 

for some constants Co £ C and po £ [0,1]. By continuity of w it is required 

2 

that lirrip^i ({]{ s 2 (l — s 2 )p~ 1 ds) = 0, which is impossible since the integrand 
is strictly positive. This proves that k = 1. 
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Finally, we establish the estimates for the semigroup. Recall that the growth 

bound wo(5), cf. [17], P- 251, for a semigroup S = (S(r)) r >o can be related 

to the spectral radius r(S(r)) of the bounded operator S(r) for each r > 0 

by the Hadamard formula. This yields wo(«S) = - logr(S(r)). From Lemma 

2 

14.91 we know that r(So(r)) < e p- 1 for all t > 0. By the Duhamel formula, 
see 02], p. 258, Prop. 2.12, 

(1 - P)S(r) = S 0 (t) + [ T So(r - T y )L , S(r , )dr' - PS(r). 

J o 

Compactness of L' and the fact that P has finite rank imply that for ev¬ 
ery r > 0 the operator (1 — P)S(t) is the sum of So(r) and a compact 
perturbation. If r((l — P)S(r)) < e~~ T for all r > 0, Eq. (14.131) follows 
immediately. If r((l — P)S(r)) > e p- 1 , then (1 — P)S(r) has a spectral 

point fj, € C with |/r| = r((l —P)S(t)) = e"~p =l+a ' >T for some a > 0. Since fx 
is not in the spectrum of So(t) it must be an eigenvalue and by the spectral 
mapping theorem for the point spectrum D3. IV.3.7, p. 277, the generator 
has an eigenvalue A with ReA = + a - I n view of the spectrum of 

L on the stable subspace, this is a contradiction if p > 3. If p < 3, then 
we know that ReA < — 1 and we infer that \p,\ < e~ T . This implies that 
r((l — P)S(t)) < e~ UJpT for all r > 0, where u> p = min{^-, 1}. This and the 
definition of the growth bound show that for every e > 0 and ji p := co p — e 
there is a constant M > 1 such that 

11(1 — P)S(r)f|| < Me~^ pT \\{l — P)f|| 

for ah f € %. The fact that A = 1 is an eigenvalue with eigenfunction g 
yields S(r)Pf = e T Pf. □ 

4.3. Nonlinear Perturbation Theory. For the rest of this section we re¬ 
strict ourselves to real valued functions. Furthermore, whenever the domain 
in the H™' 1 —norm is not indicated, it is the unit ball C R d . By B$ we 
denote the open ball in centered at the origin with radius 6 > 0. 

4.3.1. Estimates for the nonlinearity. For u = (ui,« 2 ) we define 

N(u) ('>) ; = ( N(uM) ) • 

where 

N(x) := n(c p + x) — n{c p ) — n'(c p )x, n{x ) = x\x\ p ~ , 
and c p is the constant from Eq. (11.31) . Obviously, N( 0) = IV'(0) = 0. 
Lemma 4.12. Let 5 > 0 be sufficiently small. Then 

l|N(u) — N(v)|| < (11u|| + ||v||)||u — v 
for all u, v E C Tl. 


(4.14) 
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Proof. We show that 

\\ N (u(\ ■ I)) - N(v(\ • |))||/T"d 

^ (IKI ' 1)1 \H™d + ||u(| ' \)\\H™d)\\u(\ ■ |) - «(| • |)|| H m d 


for all u, v E 0,1] that have H ™ d -norm less then 5. By density, this 
estimate can be extended to all of H™ d (M d ) and Eq. 114.141) follows. 

Note that c p > | for all p > 1. In fact, c v —>• 1 for p —>• oo and c p —>• oo as 
p —>• 1. This implies that ./V: [— —>• M is smooth. By the fundamental 

theorem of calculus, 


N(x) - N{y ) = / N'{s)ds 


= (x-y) / N'(y + s{x - y))ds, 
Jo 

for all x,y E [ — 5 , ^]- From Sobolev embedding we know that 

IMU°°(0,1) < IWI • 1)11 H m d 


(4.16) 


for all u E H md (M d ). Hence, we choose 5 > 0 so small that ||'u||l <x> (o,i) < \ 
for all ||it(| • |)||^"id < 5. Now let u, v € Cg°[0,1] satisfy this smallness 
condition. Then v(p) + s(u(p ) — v(p)) E [— 5 , for all s,p£ [0,1], The fact 
that H md (M d ) is a Banach algebra and Eq. (|4.16l) imply that 


||iV(n(|-|))-ArKM))||^ d 

< ||u(| • I) - v(| • |)||tf"*d [ || [N* O (y + s(u - v))] (I • |)|| H m d ds. 

Jo 

We estimate the integral term with Moser’s inequality, see for example [44]. 
To this end, we extend the relevant functions to the whole space. Using a 
smooth cut-off function we can construct F: R —>• M such that F is smooth, 
F = N 1 on [—4,4] and F = 0 on M \ [— |, |], The properties of N imply 
that -F(O) = 0. To extend u and v we apply Lemma IB. 21 and note that the 
extension U E C md [ 0, 00 ) of u can always be constructed in such a way that 

11^(1 ' I) llz/°°(R d ) = IMIl°°(0,1)- 

By Lemma IB.21 


11^(1 ' l)ll,ff m d(Rd) ^5 ll M (l ' 1)11 H rn d (IB^) * 

The respective extension for v is denoted by V. By Moser’s inequality, 
|| [N 1 o (y + s(u - u))] (| • |)||^ d(Bd) 

<\\[Fo(U + s(U-V))](\-\)\\ Hmdm 
<WU(\-\) + s(U(\.\)-V(\.\))\\ Hrndm 
~ I|t4(l • l)ll// m d(R d ) + ll^(| • |)||fT"d(Rd) 

~ ll u (l • 1)11 H m d (Rd) + ||v(| • | ) || jf m d (B d ) 
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for all s € [0,1]. This implies Eq. (14.151) . □ 

4.3.2. The nonlinear Cauchy problem. For u6Hwe consider integral equa¬ 
tion 

<f>(r) = S(t)u + f S(r — T , )N(^ , (r / ))dT / 

Jo 


on the Banach space 

X := {<h € C([0,oo),H) : ||$||a’ := supe MpT ||4>(r)|| < oo}. 

r>0 

Here, fi p > 0 is the constant from Lemma 14.111 In the following we denote 
by X$ the closed subspace 

X s :={<f>eX : ||$||* < <5}. 


4.3.3. Correction of the unstable behavior. We define 

r oo 

C($,u) :=Pu+ / e- T, PN(^(r'))dr', 

Jo 

and set 

K(3>, u)(r) := S(r)u + f S(r — T / )N(^'(T , ))(iT , — e r C($, u). 

Jo 

Theorem 4.13. Let 5 > 0 be sufficiently small and let c > 0 be sufficiently 
large (independent of 5). For every u € TL with ||u|| < | there exists a 
unique 4>(u) € X$ that solves the equation 

$(u) = K(4>(u), u). 

Furthermore, the map u H > 4*(u) is continuous. 

Proof. We argue along the lines of m, m- For fixed (3>,u) € X§ X Li, 
continuity of the map r i-a K(4>, u)(r): [0, oo) -A Ti follows essentially from 
the strong continuity of the semigroup. To see that K(-,u) maps X$ into 
itself for ||u|| < |, we decompose the operator according to 

K(<f>, u) = PK(4>, u) + (1 - P)K(4>, u). 

By Lemma 14.121 we have 

l|N(4>(T))||<5 2 e- 2 ^, 
for 4> £ X§ and all r > 0. Hence, 

/ OO 

e -"'||PN(4»(T'))||dT'<<5 2 e - 2 ^, 


||(l-P)K($,u)(r)|| <e-^ r ||u||+ f '')||N(T>(r / ))IMr' 

Jo 

<(S_ + S 2 )e -n P r 


and 
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Consequently, e /ipT ||K(<I>, u)(r)|| < | + 5 2 < 5 for all r > 0, given that 
c > 0 is sufficiently large and 5 > 0 is sufficiently small. For the contraction 
property of K(-,u) we use a similar decomposition and the fact that 

l|N($(r)) - N('h(r))|| < Se-^lMr) - tt(r)|| 
for any 4>, € X 5 and all r > 0 by Lemma 14.121 In particular, 

/ OO 

e r-r'(l+2 Mp ) dT / 

< 8 e- 2 ^ T \\^-^\\ x , 

and 

||(1 - P)[K(<3?,u)(r) - K(*,u)(r)]|| < Se~^ f ||<P(^) - ^(r')\\dr' 

Jo 

which implies that K(-,u) is contracting given that 5 is sufficiently small. 
An application of the Banach fixed point theorem yields the existence of 
a unique solution 4>(u) € X§. Continuity of the solution map u 1 —> 4>(u) 
follows easily from the estimate 

||K(4>, u)(r) — K(4>, v)(r)|| = ||S(r)(l — P)(u — v)|| = e _/ipT ||u — v|| 
and the fact that K(-,u) is a contraction mapping. □ 


4.3.4. The initial data operator. For R > 0 we set 

U R := H™ d x H™ d ~ l (B d R ) 


with norm defined by 


v 


I UR = IMI 


^H m d{M d R ) H^G 


I H m d ~ 




cf. (14.11) . If R = 1, then we simply use the symbol Ti , as before. In particular, 
II ' II = II ■ llw 1 - Set 


V(v,T)(p):= 

\ Tp~ 1 V2(Tp) 


, «(T):= 


( Tn ) P 1 


(£) 


P+1 
P-1 . 


and 


U(v, T) := V(v, T) + k(T) - «(T 0 ). 


Lemma 4.14. Let v £ 'H To+s for 5 > 0 sufficiently small. Then T 1 —> 
U(v, T): [To — 5, Tq + 5] — >• R is continuous. Furthermore, if ||v|+t 0 +« < <5 
then 


l|U( v , T)\\ <5 


for all T € [T 0 - 5,T 0 + 5}. 
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Proof. For simplicity we prove the result only for To = 1. The general case 
is analogous. Let v £ PL 1+S for 0 < 6 < To show continuity of the map 
T h > U(v, T) we consider the first component and estimate 

||[u(v,r)] 1 -[u(v,f)] 1 ||^ d(Bd) 

= \\Tp- 1 vi(\T ■ |) — Tp~ 1 vi(\T ■ |) + Tp~ x c p — Tp- 1 c p \\ H m d ^ B d' ) 

< ||+l(|T • |) — vi(\T • |)||H m d(B d ) + I Tp- 1 — Tp- 1 |||vi(|T • DH/p^b**) 

2 ._. 2 

-I- — TFn | 

for T, T £ [1 — 5 ,1 + (5]. Scaling implies that for all T £ [1 — 5,1 + 5] 

11^(1^' I) (B d ) ^ ll V (l ‘ I) ll_ff m d (Bj + 4 ) 

for v £ H™ d (Mf , s ). Using this and the triangle inequality we infer that for 
all vi,V! £ H™ d (Mf +s ) 

\\ v l(\T • I) — v l{\ T • I) ll/f m d (B d ) < IMI T • I) ~ Vi(\T ■ |)||H m d(B d ) 

+ ll^id^ 1 • I) — n(\ T • I)llj? m rf(B d ) + ll^id^ • I) — M\ T ■ 1)11 H m d (b^) 

~ IMI • I) ~ ^l(l • Dllir^+B^) + ll^l(l r • I) ~ Vi(\T ■ \)\\H™d(M d )- 

Since Ch+fO, 1 + <5] is dense in H™ d (Mf +s ) there is a hi £ C£°[0,1 + 5] such 
that ||ui(| • |) — hi (| • |)||^-m d ( B d ^ < e for given e > 0. From the smoothness 
of hi we infer that infer that 

hm ll^i (T-) - hi(T-)|| / / m d(B+ = 0. 

T->T 

Similar estimates can be obtained for the second component which yields 
the claimed continuity. For v £ PL 1+S , ||v||^i+$ < 5 and T £ [1 — 5,1 + 5], 

ll[U(v,T)]l||#m d ( B d) <TP- 1 \\vi{T-)\\ H rn d( ^d) + Cp\TP- 1 — l| 

~ 11^1 ll_ff m d (B^ +(5 ) + \T - 1| ~ $■ 

A similar estimate can be obtained for the second component and we infer 
that 

l|U(v,T)|| <5. 

□ 


4.3.5. Variation of the blowup time. 

Theorem 4.15. Set R := T 0 + | for c > 0 sufficiently large and 6 > 0 
sufficiently small. For every v £ % R with ||v||^h < A there exists a T £ 
[Tq — To + |] and a function <F £ X$ which satisfies 

$(r) = S(r)U(v, T) + [ S(r - r , )N($(r'))dT' (4.17) 

Jo 
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for all t > 0. In particular, $ is the unique solution of this equation in 
cm oo),H). 

Proof. Let v G U R . For 5 and c chosen appropriately, the smallness con¬ 
dition for v and Lemma 14.141 imply that U(v,T) satisfies the assumptions 
of Theorem 14.131 for all T G [To — |, To + |]. Hence, for every such T there 
exists a <h;r := <&(U(v, T)) G Xg satisfying 

$ T (r) = S(r)U(v,T) + [ S(r - T , )N($ T (r'))dT / 

Jo 

-e r C($ T ,U (v,T)) 


for all t > 0. We show that there is a T v G [To — |,To + |] such that 
C(<f>T v ,U(v,T v )) = 0. Note that rg C = (g), where g is the symmetry 
mode from Lemma 14.101 Hence, it suffices to show that 


(c($ Tv ,u(v,r v ))|g) =o. 


We find that 

roo 


poo 

/ e" T PN($ T (r))d7 

Jo 


g J ^ l|g 


< 


poo 

/ e- r ||PN(<i> T (T))||d7 

Jo 

poo 

Jo 


(rWdT<5 z . 


The key observation is that 


<9tk(T)|t=To = a g 


for a constant a > 0. By Taylor expansion of n we get 

U(v, T) = V(v, T) + (T - T 0 )ag + (T - T 0 ) 2 q T 
for all T G [To — |,To + |] and for some remainder term q^. Thus, 

(PU(v, T)|g) = (PV(v, T)|g) + a(T - T 0 )||g|| 2 + O(g) 

= a(T-T 0 )\\g\\ 2 + O(j I ) 

by definition of V and the smallness of v in U To+ c. We notice that the 
order terms depend continuously on T. Summing up, we obtain that the 
equation 

(C($T, U(v, T))|g) = a(T - T 0 )||g|| 2 + O(^) + 0(5 2 ) = 0 


is equivalent to 


T = T 0 + F(T) 


where F(T) = O(-§?)+0(5 2 ). For c sufficiently large and 6 = 5(c) sufficiently 
small we get |T(T)| < Hence, the continuous function T i —> Tq + F(T) 
maps the interval [To — To + |] to itself and has thus a fixed point at some 
T = T v . We therefore obtain a solution <f>r v € Xg of the original equation 
(I4.17p . For the uniqueness of the solution in C([0, oo),'H) we refer the reader 
to the proof of Theorem 4.11 in [10] . □ 
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4.3.6. Proof of Theorem \1.4\ Choose 5, c > 0 such that Theorem 14. 151 holds 
and set S' = S/c. Let (uo,rti) € H^~ x H^~(M. d ) be radial functions, i.e., 


Oo^i) = 

For (f,g) 


(&o(l ' |),«i(| • I)), that satisfy 
||Oo,'«l) - U To 


< 2+1 < 2-1 


< S -. 

— c 


= (iio,ui) — ut 0 [0], cf. Eq. (|3.2I) . this assumption implies that 


\\(fi9 )\\h t o+ s / c ^ j?- 


Hence, v := (f,g) satisfies the assumptions of Theorem 14.151 We infer that 
there exists a T E [To — 5, To + <5] such that Eq. (14.171) has unique solution 
<f> E C([0, oo),72) with 


||$(r)|| < 5e~^ T , Vr > 0. 


Hence, T = <f> + c p is a solution of Eq. (13.31) (in the Duhamel sense) with 
initial data T(0) = U(v,T ) + c p . Consequently, 

u(t,x) = (T — t) ^ rT '0i(— log(T — t) +logr, J^L) 
is a radial solution of the original wave equation (11.11) with initial data 

u( 0,x) = i(0, = ut 0 (0, x) + f(\x\) = u 0 (x) 

_ p+1 | | 

d t u(0,x) = +) = d t u To ( 0,x) + g(\x |) = ui(x) 

for x € By. Furthermore, u satisfies the estimates 
( T—t)p _1 ^ \\u(t, •) — ux{t, •)_ t ) 

= (T - tyi +k ||^r(- log(T - t) + logT, - c p \\ kk{M d_ t) 

= IIVh(- log(T - 2) + log T, | • |) - c p \\ tik{Bd) 

= ||pi(-log(T-f) + logT, | • |)||jjfc (B d) 

~ l|Pi( — l°g(T — t) + logT, | • |)||/f-™d(B d ) 

<||4>(-log(T-t) + logT)||<(T-t)^. 


for k = 0,... , The bounds for the time derivative of the solution follow 
accordingly. 


5. Improvement of the topology - Proof of Theorem 11.51 

The verification of Theorem 11.51 is analogous to the proof of Theorem 11.41 
and we only discuss the main arguments. With the definitions of Section 
IQ we introduce the product space % := H™ d ~ 1 X H™ d 2 (B d ) with norm 

ll U lll — ll^ld ' I ) 11 — 1 (IB rf ) "F 11^2(1 ' \)\\f[m d -2 ^dy 
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5.0.7. Time evolution for the linearized problem. We proceed as in Section 
HEED Since most proofs are similar or can even be copied verbatim we only 
sketch the main steps and point out differences. With the same notation as 
in Section IP we define 

\\M\d := ll^' Ul llj/i(0,l) + 11IIL 2 (0,1)- 

Lemma 5.1. We have 

lllulll 2 =* \\D d Ul\\ 2 H i m + 11^21112(0^) ^ W\u\\\ 2 d 

for all u € 77. Furthermore, D d u\ £ C[ 0,1] and D d ui(0) = 0. 

Proof. The equivalence of the parts involving only rq follows from Lemma 
BZD and the proof of Lemma 14.21 For the second component the same meth¬ 
ods can be used to show that 

m d -2 

INI-DllWw^ £ IK')” +1 “ ( ” , llw(0,i) = 11^11^(0,1) 

72—0 

for all u £ H™ d ~ 2 (M d ). The properties of D d u\ are a consequence of the 
density of C^°[0,1] in PL and the embedding H l ( 0,1) ^ L°°(0,1). □ 

We note that Eq. (14.31) and Eq. (14.41) of Lemma ITT 1 hold for all u € H md ^ 2 n 
C md ~ 1(0,1) and u £ FT"'* -1 n C md { 0,1), respectively. For the moment, we 
leave the value of p unspecified and consider Lo as defined in Eq. ()4.7|) on 
the domain 

D(L 0 ) := {u £ PL n C°°(0, l) 2 : D d u 1 £ C 2 [0, 1], 

DdU2 £ C ,1 [0,1], ( DdU2)(0) = 0} . 

The perturbation L' is defined as in Eq. (14.101) and Lemma 14.81 holds with 
PL replaced by PL. It is easy to check that 

Re(D d {L 0 u)i\D d u 1 ) ^ 1(0 1} 

+ Re(D d (L 0 u)2\D d u 2 ) L2{0A) < (5 - ^= T )|||u|||| > 

for all u £ D(Lq) and that rg(yu — Lo) is dense in PL for /jl = 1 — . With 

L 0 denoting the closure of Lo, we use the same arguments as in Section 14.2.11 
to infer that L := Lo + L', D(L) = D(Lo), generates a strongly-continuous 
semigroup (S(r)) T> o of bounded operators on PL. As an analogue to Lemma 
14.101 we obtain the following result for the spectrum of L, where we fixp = 3. 

Lemma 5.2. Let p = 3. If X € <r(L) and ReA > then A = 1. Further¬ 
more, it is an eigenvalue and ker(l — L) = (g), where g = (1,2). 

Proof. The assumptions on A imply that A fL ct(Lq) and the fact that A is an 
eigenvalue follows from the compactness of Lb By definition, if u £ PL, then 
u\ and U 2 are 1-times, respectively, m d — 2-times weakly differentiable. 
Sobolev embedding yields ui € C md ~ 2 [5, 1], U 2 € C md ^[5.T\ for arbitrary 
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5 > 0. For d > 7, this already implies that eigenfunctions corresponding 
to the eigenvalue A satisfy the equation (A — L)u = 0 in a classical sense 
on (0,1). For d = 5, one can use the definition of the closure to check that 
Ui £ H^ oc ( 0,1) n C 2 (0,1), U 2 € H 2 oc ( 0,1) n C' 1 (0,1) if u S D(L). Hence, 
for all d > 5 odd, the first component of an eigenfunction solves Eq. f)4.1 1 [i 
in a classical sense on (0,1). By smoothness of the coefficients on the open 
interval, we get that u\ £ C 00 ^, 1) and the application of Dd shows that 
w := DdU\ solves Eq. (14.1211 on (0,1). By Lemma l5Tl w £ H 1 ( 0,1) n C[0,1] 
and w(0) = 0. We now argue as in the proof of Lemma Pi. 101 cf. also pH], to 
infer that A = 1 and that the corresponding eigenspace is spanned by g. □ 

With similar arguments as in the proof of Lemma 14.111 we obtain the fol¬ 
lowing result. 

Lemma 5.3. Let p = 3. There exists a projection P C B(TL), rgP = (g), 
that commutes with S(r) and 

S(r)Pf = e r Pf 

for all f £ TL and all t > 0. Moreover, 

|||(1 - P)S(r)f||| < Me“ ( 5- £ ) T |||(1 - P)f||| 

for all f £ 74, r > 0, some M > 1 and some small s > 0. 

5.0.8. Lipschitz estimates for the nonlinearity. In the following, B denotes 
the unit ball in TL. For p = 3, the nonlinear remainder is given by 

N <"> ;= ( a? + 3^; ) ■ 

Lemma 5.4. The operator N: TL -A TL satisfies 

|||N(u) — N(v)HI < (HI u||| + HIv|||)HI u — v||| 
for all u, v £ B C TL. 

Proof. For u £ M d ), rrid — 1 = we set u(£j) := u(|£|), £ £ M d . In 

the following we do not indicate the domain in the Sobolev norms, since it 
is always the unit ball M d C M d . The Sobolev embedding W q+m ’ 2 ^A W q,q 
for 2 < q < d 2 ^ m implies that 

\\d a u\\Li < \\u\\ H m d -i 

for a £ N d , 0 < |a| < and 2 < q < We first consider the cubic 

part of the nonlinearity. To estimate the L 2 -part we use Holder’s inequality 
with qi = q 2 = 2d, ± ^ = |, to show that 

||S 3 — v 3 \\lz = IK'S — v){u 2 + v 2 + uv )\\ L 2 

4 11U — S 11 L®1 ( 11S 11 fq2 + 11 V 11 \q2 + 11 S 11 £92 11 V \ \ LI 2 ) 
iS (IH|^m d -l + ||h||^m d -l)||'U — '0||^m d -l. 
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For higher order derivatives we have to estimate terms of the form 

dP(u — v )d a ^ l3 (u 2 + v 2 + uv ), 

for 0 < /3 < a. For |a| = J-^-, ft = a, we apply again Holder’s inequality 
and Sobolev embedding to get for example 

\\(d a u - d a v)ii 1 \\ L 2 < \\d a (u — v)\\l«i IHIi? 2 < \\u ~ v\\ H m d -i\\u\\ 2 H m d -i 

for qi = q 2 = 2 d. Since d a u 2 is equal to a sum of terms of the form 
d ai ud a2 u, where oq + a 2 = a, we infer that for (3 = 0, 

||(u - v)d ai ud a2 u \\ L 2 < ||u - v\\ L n \\d ai u\\ L i 2 \\d a2 u\\ Lq3 

< HU — '0||^m d - 1 ||'u||^ md _ 1 , 

where qq = 2d, q 2 = Q 3 = i + lf Q2 | ^ E/=i ^ All other terms can 

be estimated similarly. 

For the quadratic part of the nonlinearity we set for example qq = 2d, q 2 = 
/FT ’ to g et 

||(w - v)(u + v )\\ L 2 < ||u - v\\ L n II u + v\\ L Q 2 

< \\U + v\\ H m d -l\\u — v\\ H m d -l 

or, for |a| = qi = 2d, q 2 = 

\\(u - v)d a u \\ L 2 < ||ft - v\\ L n \\d a u\\ L n < \\u + v\\ H m d -i\\u\\ H m d -i. 

Estimates for the remaining terms follow from similar considerations. □ 

With Lemma 15.31 and Lemma 15.41 Theorem 11.51 follows by proceeding as 
above, starting with Section 14.3.21 

Appendix A. Hardy’s inequality 

Lemma A.l. Let a € N. Assume that f € 0°°[0,1] satisfies 0) = 0 
for j = 0,..., a — 1. Then, 

ll(-)"”/IL> < IIC-)-“ +1 /'lli=- 

Proof. For / = 0, the assertion is trivial. Let / 0. We use integration 

by parts, l’Hospital’s rule and the Cauchy-Schwarz inequality to obtain the 
estimate 

f 1 p- 2a \f(p)\ 2 dp< lim (|2a- l|-V 2a+1 |/(p)| 2 ) 

Jo o 

+ |2a — lp 1 [ p~ 2a+1 {f'{p)f{p) + f{p)f'{p))dp 

Jo 

< f 1 p- 2a+1 Mf'(p)7(p)]d(< f P - M \f'mW)\dp 

Jo Jo 

< [j\- 2a \f{p)\ 2 dp^ 1 (f\- 2a+2 \f'(p)\ 2 dp) 


1/2 
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This implies the claim. 

Appendix B. Proof of Lemma 14.11 


□ 


It suffices to show that the claimed inequality holds for all u € C e °°[0,l] 2 . 
By density this can be extended to all of LL. In the following we set 


X7 n •— 


A 


n/2 


for n even, 


rad ' | d_ A {n-l)/2 fornoddi 


where A P u(p) = u"(p) + ^-J-u'(p). To abbreviate the notation we define 


m d 


|u|| li IMll 2 ro.D + E IK')” lu(n) lli 2 (o,i)’ 

m d -1 


Mls 2 


n— 1 

r.i n 7/.( n )ii 2 


(B.l) 


£ ll (■)'■«' - iii» ( o,i)- 


n =0 


Lemma B.l. Let u € C%°[0, 1]. Then 

IMIs, < IHI • 1)11 and IMIs 2 ^ IKI • DIIh"**- 1 

Proof. We prove the first estimate. Let V := (d\, ..., dd) T ■ Then 


m d r i 

£ 

n=0 


m d 


E/ P d_1 |V?ad«(p)| 2 dp^E / l Vr MI£l)| 2 ^ 


< E -1)11^0-) = HI 

\a\<m d 

In view of this, it suffices to show that 

m d r l 


I 2 H m d 


\u\ 


2 < 
El 


" v a /*! 

E / ^ d_1 | V ?ad“(p)| 2d ^ 

n=0 J ° 


(B.2) 


First, observe that 

m d 

imi1 2 (o,i)+E lu(n) ni 2 (o,i) 


n=1 


rl m d -l I 

II' u IIl 2 (o,i) + / P d ~ 1 \V™* i u(p)\ 2d p+ E / p 2(n " 1 V n) (,9)| 2 dp 

Jo n=1 Jo 

rl rl 171 d 1 

/ Kp) | 2 + / p d_1 |V^li(p)| 2 dp+ E / p 2(n_1) l v ! 1 ad^(p)l 2 rfp- 

JO JO , JO 


< 


< 

rsj 


71=1 “'° 


We show that 


-i 


"id /■! 


P 2 (n_ 1 ) |V r ^(p)I^P^E / A 1 | VJ rad«(P)| 2 ^ (B. 3 ) 

1=0 ^ 
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for all n = 1,..., — 1. Recall that for radial functions the trace theorem, 

cf. for example m, p. 258, implies that 


|V? ad R. 


( 1 )| 2 < f 1 P d - l \K^(p)\ 2 dp+ C p d - l \V^u{p)\ 2 dp 
Jo Jo 


for all n = 0,... — 1. Assume is odd and let n = — 1. Then 

integration by parts and the Cauchy inequality imply that 
-1 


L 


< 

rsj 


< 


p d ~ 3 \K^Mp)\ 2 d P 

ra _ 1 «Ki) 2 + f 1 p d - 2 m~Mpw?>i(p)\dp 

Jo 

+ I I' p d ~ l \KMp)\ 2 dp + e / 1 /- 3 |V^- 1 u( /0 )| 2 dp, 
e Jo Jo 

for any e > 0. For mu even one can easily check that the function (■) <i_1 V™j _1 r 
satisfies the assumption of Lemma I A. 1 1 for a = (d + l)/2. Hence, Hardy’s 
inequality can be applied to obtain 
/•l 


p d ~ 3 


< 


l u{p)\ 2 dp = f 1 p- d - 1 |/- 1 [V r ^- 2 ^] / (p)| 2 ^ 
Jo 


d+l 


,rf— 1 r X7 m d~ 2 n,V 


L [v 


rad 


U 


dp 


< 


p d ~ l \KMp)\ Z dp. 


Now these arguments can be iterated to get Eq. (IB. 31) . Finally, 

[ \u{p)\ 2 dp < \u{l)\ 2 + f p\u(p)\\u'(p)\dp 
Jo Jo 


< 


«(l)l 2 + 


if 

e Jo 


\u'{p)\ 2 dp + £ / \u(p)\ 2 dp, 


for any e > 0 and the first line of Eq. (IB. 2D follows. The second estimate in 
Lemma lB.il can be obtained analogously. □ 


Lemma B.2. Let u G Cg°[0,1]. There exists a compactly supported function 
U G C md [0, oo) such that U(p ) = u(p) for all p G [0,1] and 

ll“(l ' I) (B d ) < 11^(1 ' I) lllL m d (R d ) ~ Iklbi- 

Similarly, one can construct an extension U G C md ~ l [ 0, oo) such that U(p) = 
u{p) on [0,1] and 

11^(1 ’ l)lllf m d _1 (B d ) ~ 11^(1 ' l)ll// m d- 1 (R d ) ^ IM|e 2 - 

Proof. For m G N, let / G C°°[0,1] and let <p G C 00 ^, oo) be a monotonically 
decreasing function such that ip = 1 on [0, |] and ip = 0 on [|, oo). We define 

/(p) for p G [0,1] 

V(p)hm(p) for pG (1,3/2) 

0 for p G [3/2, oo) 


£mf(p) ■= < 
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where 


h m {p) := < 


(m ^ 1)/2 2(n-1) 2n 

-f(2~p)+ E (2n)' /(2rt)(1) misodd ’ 


n =0 


m/2 Of n _ 1 ^2n—1 

A2-rt + E Ljn, ™ is 


even. 


^ (2n — 1)! 

i=i v ' 

Then, £ m f G C m [0,oo) and £ m f |[o,i] = /• We define U := £ md u- Our aim 
is to prove the estimate 

||t/(| • Dlll^j + ||V"*tf(| • |)||i 2(Rd) < Ml,. (B.4) 

Given that Eq. (IB.4I) holds, we can use the fact that H m (M. d ) can be equiv¬ 
alently defined in terms of the Fourier transform to infer that 

hi • m 2 H m d{Bd) < m ■ dii^,^) < \\(-r d nu(\ • i)]iii 2(Rd) 

£ 11^(1 • DIIlw + l|v md c/(| • Dllia^) < IMIlv 

If md is odd we estimate 


roc 

l|C(| ' Dlliijna) < J p d -'lU(ptdp 

= Jo pd ~ 1 ^ p ^ dp + / 


(m d -l)/2 


< 


< 


Il 2 (o,i) 


+ 


l-[o,oo)( I* P d 1 H2-p)\ 2 dp+ E l« (2n) (l)| 2 ) 

' Jl n =0 ' 


(m d -l)/2 


Il 2 (o,i) + / \ u (p)\ 2 dp+ E/ k (2n ^(l)l" 


< 


\u\ 


Si> 


72—0 


where we used the fact that 




i“ w (i)i 2 s s ii-iihp,!) + E ii(')' _ 1 “ 0 ) ni.(o,i) 

3 = 1 


for n = 0,..., — 1. For the derivative we get 

roc r 1 

/ /- 1 |v^c/ 1 (p)| 2 dp= / p^-'Iv^Xp)! 2 ^ 

JO JO 

+ / 2 Pd_1 l V r™d[^)WP)]| 2 ^- 

To bound the first integral we exploit the fact that for d > 5 odd and m G N 
there exist constants c\ L ' m ' 1 g R such that 


V£d -(rt = E c »'”V~”‘« ( ” ) (p)- (n.5) 

72 — 1 
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Hence, 


r 1 rl 

Jo n=1 J 0 

It remains to estimate the second term. It follows from Eq. (IB.5D and the 
Leibniz rule that there exist constants G R such that 

,L iJ 

m d n 

V”' lv(p)h mj (p )] “EE'SA'fMWieif). 

n= 1j=0 


Hence, 


P d 1 TO [^(P)^ d (p)]| 2 ^ 


m d n „ 3 

SEE f 


„2n—2 


n=lj=0 
/ m d 


'1 


^ n - j \p)hiil(p) dp 


< 

r^j 


EM”’" ? 


L°°[0,oo) 


^n=0 


' m d 

E 

^n=0 ' 


„d-l 


1 




< V f 

~ / / 


u 


n=0 


'1 


(m d —1)/2 

( n )(2- /0 )| 2 d / 9+ M 2 ")(l)| 2 

71—0 


m d r i 

/ l^ (n) (p)| 2rf P + II«IIe 1 ^ IMllv 

^—n J tt 


< 


n=0 2 


If is even the proof works similarly. In fact, the extension was constructed 
in such a way, that the boundary terms involve only derivatives that can be 
bounded by the E-norm. The proof for the second estimate is analogous. □ 


Appendix C. Proof of Lemma 14721 

Again, it suffices to prove the inequality for all u G Cf°[0, l] 2 . We split 
the proof into several lemmas and use the result of Lemma 14.11 With the 
definition (IB. ID we set 


||u|| E := ll'UlHsi + ||U2 ||s 2 

for u = (iii, U 2 ) G C“[0, l] 2 . 

Lemma C.l. We have that ||u||_d < ||u||s, for all u G C£°[0, l] 2 . 
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Proof. With Eq. (14.21) and the triangle inequality we immediately obtain 

m d 


IIAKI 


H 2 ( 0,1) ~ I 
n=l 


j-l ,(n )||2 


n 


lll 2 (0,l)’ 


\\ D d u 2 1| i4i( 0 q) 


m d -1 

< IIAW„,( n )l|2 


~ E iko 

n=0 




IIl 2 (0,1)- 


We use the Sobolev embedding i^ 1 (0,1) L°°(0,1) and the fact that 

D d U 2 (0) = 0 to infer that 

IKK'U) + D d u 2 ( 1)| 2 < \[D dUl }'(l)\ 2 + \D d u 2 (l)\ 2 

~ ll[^d M l] IIiA 1 (0,1) + llAKI^o,!)' 

Now, 

rn^-l 

II IKK llz,2 ( 0jl ) ~ II ^ 11 Zy 2 (0,1) 


n =0 




~ II U iHl 2 (0,1) + ^ 


sre-1 (n)ii2 

1 u i IIl 2 (o,i)’ 


n=l 


□ 


which implies the claim. 

Lemma C.2. For all u € Cf?\ 0, l] 2 the following inequality holds 

11 (0,1) HKKI/f 1 ) 0 .!) ^5 II u IId- 

Proof. Set Wj := D d Uj for j = 1,2. Since Wj( 0) = 0 for j = 1,2, the above 
inequality is true if 

lKlli 2 (0,l) < K(l) + ^2(l)| 2 + IK 11^2(0,1) + ||«4|lia (0>1) . 

By the fundamental theorem of calculus fj wf(s)ds = w[( 1) — w\ (p) and 
Ip w' 2 {s)ds = 102 ( 1 ) — W 2 (p)- Hence, 

\w[{p) + w 2 (p)\ < K(i) + 102 ( 1)1 +1 IK 11L 2 ( 0 , 1 ) + IIKIIl 2 (o,i)- 

Using this together with Sobolev embedding we obtain 

K(P)I < IK(p)+^ 2 (p)| + KKI 

< |K(i) + ^ 2 (i)| + IIKIU 2 + IIKIIl 2 (o,i)- 

Squaring and integrating implies the claim. □ 

Lemma l4. 21 follows from Lemma f4. II in combination with the following result. 
Lemma C.3. We have that 


for all u € 1] 


1 2 < 
Is ^ 


IK n llK(0,l) + IKKIffl(0,l)> 
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Proof. We show that 

ll(-)"« < " ) IR=(o,i) < (c.i) 

for n = 0,..., nid — 1 and all u £ Cg°[0,1], and 

||(.)"- 1 « ( ” ) || 1 .(o, 1) <C||B J »|| H3( o, 1) (C.2) 

for n = 1 by using the fact that KdDd = I on C“[0,1]. Let 

u € C£°[0,1] and set w := DdU. Then, w £ C°°[0,1], w^ 2n \0) = 0, n £ No, 
and KdW = u. One can easily check that there exist constants a n j,a n j £ R 
such that 

n (d—3)/2 

p n {K d w) ( - n \p) = E a n .jK d _ 2 ]W = y a n pp~ 2:, ~ l Vw(p) 

3 =0 j=(d—3)/2—n 

for n = 0,..., nid — 2 = (d — 3)/2. Furthermore, 

{d- 3)/2 

p rn d- 1 (K d w)( md ~ 1 \p) = 'y aijp- 2 ^- 1 ^ w(p) + w' {p). 

3=0 


Since w( 0) = 0 and 

(o) = o 

for n £ Nq, k = 0,... , 2n, repeated application of Hardy’s inequality yields 


f 


—An—2 1 


|/C n u;(p)| 2 (ip < C n / |u;'( / 9)| 2 (i / 9 


for n € No and some constante C n > 0. This in particular implies that 

/V(W n) (p)| 2 ^< Clw'iptfdp 

Jo Jo 

which proves Eq. (1C.II) . By the fundamental theorem of calculus 

n s 

w"(t)dtds. 

Upon setting Vw(p) := w(s)ds we infer that 

K d w{p ) = I< d V 2 w"(p ) + k d w'(0 ) 


for some constant kd > 0. Using this we obtain for n = 1,..., nid — 2 and 
constants /3 n ,/3 n £ R, 

n 

p n ~ i {Kdw) ( ' n) (p) = y^p-'K d _ 2 y w "{ P ) 

3=0 

(d- 3)/2 

= E Pn,jp- 2j - 2 1C j V 2 w"(p). 

j=(d—3)/2—n 






Furthermore, 


(d— 3)/2 

^- 2 (i(» K - 1 )(/))= 7i p- 2j '- 2 /C' J V 2 '«;"(p) + ^ 1 V^"(p) 


p m d~ l (K d w)( md \p) = ^ j jp - 2j - 2 JC j V 2 w"(p) 

3=0 

+ 7 dP~ 1 Vw"(p) + w"(p) 

for 7 ?, 7 j, € M. By repeated application of Hardy’s inequality we can now 

show that 

[ p~ 4n ~ 4 \)C n V 2 w"(p))\ 2 dp < C n [ \w"(p)\ 2 dp 
Jo Jo 

for n € N and a constant C n > 0. This implies Eq. ()C.2I) . □ 
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